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Abstract 

An analytical approach to network dynamics is used to show that 
when agents copy their state randomly the network arrives to a station- 
ary regime in which the distribution of states is independent of the degree. 
The effects of network topology on the process are characterized introduc- 
ing a quantity called influence and studying its behavior for scale-free and 
random networks. We show that for this model degree averaged quan- 
tities are constant in time regardless of the number of states involved. 
Keywords: Complex Networks, State Dynamics, SocioPhysics, Social Dy- 
namics, Opinion Formation. 

1 Introduction 

Complex network theory has nourished as an effort to explain the main charac- 
teristics of discrete interacting systems. The first studies were able to describe 
and explain their topology by showing that different thermodynamical quantities 
are needed to characterize a particular network, such as its degree distribution 
P 1213 EE], community structures H3 0, resilience |S] and motifs [5]. Re- 
cent efforts have also studied quantities spreading on complex networks, such as 
computer viruses |1U1 111! 1121 1131 in technological networks and information 
in their social counterpart |15) . 

On the other hand, state dynamics has been studied disguised as opinion 
formation processes, which have been modelled using approaches including cel- 
lular automata |16) . spin chains |17l I18| and opinion drifts over a continuous 
opinion space t 19: , without incorporating topology. Recently, network topology 
has become more relevant and computer simulations have been performed in 
this direction |2(J1 1211 12"2] , though a systematic way to evaluate the topological 
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effects on the dynamical processes taking place on networks is still missing. 

In this paper we study the effects of topology on a copy mechanism taking 
place over a network with a non-trivial topology and an arbitrary number of 
states. To achieve this, we group agents with similar topological characteristics 
and study the density of the evolving quantities over these groups by considering 
the interactions among them. 

We will first introduce the model together with its associated notation and then 
solve it analytically. We will then study how network topology affects the sim- 
ple dynamical process introduced below and show that the system arrives to a 
configuration in which the distribution of states is independent of the degree. 

2 The Model 

We consider a process in which at each time step an agent changes its state by 
copying one randomly from one of its immediate neighbors. The probability 
that this occurs is the product between the probability that the agent is chosen, 
times the fraction of neighbors in a given state, including the agent, that are ac- 
quainted by it. We will consider a directed network with an adjustable topology 
as the substrate in which these dynamical processes will occur, and attempt to 
characterize the influence of topology upon them. We assume that changes in 
the network topology take place at larger timescales than the relaxation times 
of the system. 



2.1 Definitions 

We begin by grouping agents according to their degree into aggregates we call 
guilds. This reduces the problem to one as big as the number of guilds used to 
calculate the approximation, this depends on the desired degree resolution, but 
it is always considerable smaller than the number of agents. Agents in a specific 
guild can be in any of the available states. We introduce the density pi a as the 
number of agents in the i th guild which pertain to the a state, over the number 
of agents in the network. We define guild size pi — ^ Q Pi a as the fraction of 
agents constituting that guild and fractional abundance as the fraction of agents 
on a particular state in the entire network p a — '^2 i Pia- From now on, latin 
indices will refer to guilds, while greek indices will label states. 

We assume that the linking patterns of agents in the same guild are similar. 
In other words, we will model a network in which agents from the same guild 
are linked to a similar number of agents in other guilds, although not necessarily 
the same. This motivates us to introduce Pij as the probability that an agent 
in the i th guild is linked to an agent in the j th guild. We will refer to these 
quantities as linking probabilities and say that someone is linked to someone else 
if the linking agent has a personal or media based knowledge of the linked one. 
An example is shown in figure (Q. Three equal size guilds are clearly distin- 
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Figure 1: Undirected network with eleven nodes in which 3 guilds are distin- 
guishable. Nodes carry guilds as its labels, the first guild has 5 agents with a 
degree equal to 2. The second guild is formed by 2 agents with a degree equal 
to 4. The third guild is formed by 4 agents with a degree equal to 1. 



guishable and acquaintanceship probabilities are given by 

/ 2/5 1/5 
P= 1/5 1 1/2 
V 1/2 

2.2 Topologically Considerate Mean-Field Approach 

In the large network limit the rate of change for a state in a given guild is equal 
to the difference between the probabilities of wining and losing an agent. The 
probability that an agent in the i th guild turns into the a th state is 



n£ = - ^ , (i) 



while the probability that an agent in the i th guild and the a th state changes 
its state to a diffrenet one is 

Pia PijPjv 

nr Q - • (2) 



We simplify the notation introducing T as the matrix representing the ex- 
pected fraction of agents in a particular state known by an agent in a particular 
guild on the entire network, 

^ la ^ ^ Pj j Pja • (3) 
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Combining eqns. QJ, (J2J and we can write the temporal variation of 
state abundance for a given guild as 

where = Ti M represents the fraction of agents acquainted by the i th 
guild 1 . The extension of the restricted sums in equations and J3J to include 
all states has no effect on the system dynamics because the additional terms 
cancel out when we place them on eqn. Time units represent the network's 
natural update time, which is the number of time steps equal to the population 
size in a random or sequential updating processes. 



2.3 Analytical Solution 

An analytical solution can be obtained by re-writing eqn. 10} and performing 
the substitution pi a — Uiae^ 1 . This transforms eqn. into an eigenvalue 
equation for the vector U a — {U± a , TJ% a ^ ■ • •). The problem then reduces to 
finding the solution of 

d t U a = MU a , (5) 

where My- = Y~Pij. From eqn. JSJ) it follows that the eigenfunctions of M 
depend exponentially on time with a decay rate A, where A is the associated 
eigenvalue, this agrees with the results presented in j2SH2E|- A general solution 
can be written in terms of this basis as 

V 

in which V? is the i th element of the rf h eigenvector of the matrix M and A ?) is 
its associated eigenvalue. It is important to notice that M does not depend on 
the state, implying that the state dependence is introduced through the initial 
conditions only, represented by C? a . 



2.3.1 System Equilibria 

The steady states are characterized by the eigenvector associated with the eigen- 
value A = 1 . By direct calculation one can show that M always has an eigenvalue 
equal to 1 associated with guild sizes as its eigenvector. One has 

Ew \ ~" PiPijPj \ •> PijPj 

m *oPj = E —f— = p*1^ ~y~ = pi - 

3 3 3 1 

1 From now on we will assume that quantities with a missing index have been added over 
it, for example: p a = ^ pi a 
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We have observed numerically that the remaining eigenvalues have a real part 
smaller than 1 and lie inside the unit circle of the complex plane, a formal 
demonstration of this remains open. 

The mean-field approach presented above was compared with a stochastic 
simulation on a 1000 agents network in which two equal size guilds were distin- 
guishable (fig. 121 )• The acquaintanceship probabilities were chosen such that 
all agents acquainted everyone in their respective guild, while the ones in the 
first guild acquainted an average of 3/5 of the agents in the second guild. On 
the contrary, agents in the second guild acquainted 1/5 of the agents in the first 
one. The simulation shows how accurate the mean field model is and reveals 
that the standard deviation tends to stabilize as the system relaxes to equilib- 
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Figure 2: The continuous line represents the mean- field solution of a two guild 
network with 1000 agents, while circles represent an average of 50 stochastic 
simulations. Error bars indicate one standard deviation. In each simulation the 
network was chosen randomly constrained to satisfy the conditions presented in 
the text. 



3 Degree Distribution 

We define the in degree of an agent pertaining to a specific guild as the fraction 
of agents in the total sample that it knows. It, therefore represents the fraction 
of the network exerting a direct influence on it. On the other hand, an agent's 
out degree is the fraction of the entire network that knows the agent, these are 
the ones directly influenced by it. This definition of degree implies that peer 
influence flows opposite to acquaintanceship. According to this, the in and out 
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degree's in our model are represented by 



thus degree distribution can be imposed by considering that the fraction of 
agents with the described degree is equal to guild size pi. We can approach a 
complex network with an out degree distribution satisfying a function / by the 
solutions of 



Pi = f[Nj2 P iiPi) (8) 

j 

In the case of scale-free networks eqn. © becomes 

I>,W) 7 > ( g ) 



where constant factors like normalization or network size have been absorbed 
by the elements of P, noting that equation is zero degree homogeneous in 
P. For the case of a random network, degree distribution has an exponential 
tail that can be approximated by the solutions of 



-1lYl P jiPj) 



Pi = A exp -vY] PjiPj (10) 



as parameters. Here rj characterizes the rate of exponential decay and A is a 
normalization factor. 



4 Guild Influence 

4.1 Scale-Free Networks 

To capture the influence exerted by the topology on the dynamical process, we 
define guild influence U as the equilibrium fractional abundance that a state 
will achieve assuming that it was initially occupying, only and entirely, the i th 
guild (p a (t = 0) = p^, 

h = Um Em®. 

t^QO pi 

This quantity represents the average number of copies that an agent of a 
given guild creates during the entire process, or an effective infection rate proper 
of the guild. This is a formalization of what Wu et. al.[2Sj and Sucheki et. al.|25| 
have done. Influence measures the affect of the group of agents with a certain 
degree, generalizing the bias depicted in the two state systems studied by them. 

We studied guild influence on a scale-free network by choosing guild sizes 
and finding acquaintanceships probabilities for a given scale-free exponent ac- 
cording to eqn. Figure JSK) shows that the influence of the highest degree 
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Guild Clustering log a 

Figure 3: (A) Guild influence against scale- free exponent 7 for a 5 guild network 
calculated through numerical integration of Eq. J3J. From top to bottom, 
guild sizes were chosen as 9 x 10 -1 , 9 x 10~ 2 , 9 x 10~ 3 , 9 x 10~ 4 , and 9 x 
10~ 5 . (B)Influence structure of a random network against the inverse of its 
characteristic length 77, for a 10 guild network. From top to bottom guild sizes 
were chosen as 0.026, 0.032, 0.041, 0.053, 0.068, 0.088, 0.114, 0.147, 0.189 and 
0.245. (C) Influence for a ten guild network in which each guild had a particular 
clustering coefficient. Different lines represent linearly spaced variations of the 
constant a in the 0.01-0.09 interval from lowest to highest slope respectively. 
The constant b was taken as 0.1 and 10 guilds were used. (D) Characteristic 
times against inter-guild connectivity for values of b chosen linearly spaced from 
top to bottom in the 0.1-0.01 interval 



guild increases with decreasing 7, until it reaches a steady value which is equal 
to the ratio between the system size and the number of agents in the highest 
degree guild. This is a consequence of the increasing fame of the highest degree 
guild that occurs for small values of 7. In the large exponent limit the influence 
of all guilds approach unity, meaning that there are no high out degree agents 
exerting an important influence on the system. 



4.1.1 Exponential Networks 

An example involving a exponential network approximated by 10 guilds is shown 
in fig. l|33), the procedure used to determine the parameters being analogous to 
the one performed for the scale-free network. The figure shows that in a random 
network high degree guilds exert a mild influence compared to that occurring 
on its scale-free counterpart. This is a consequence of the reduced population 
of high degree guilds. 
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4.2 Clustering 



The topological influence of clustering was studied by choosing acquaintanceship 
probabilities as 

Pij = a + Sij(b ■ i — a) with b > a 

and equal guild sizes for all the network, where 5ij is the Kronecker delta. This 
model represents a random network in which agents of the i th guild have a prob- 
ability a of being acquainted to agents in other guilds. As usual, we shall call 
this connectivity. On the contrary, the probability b ■ i to acquaint someone from 
its own guild together with the condition b > a defines clustered groups with 
an internal connectivity that increases linearly throughout the network. This is 
similar to communities in the sense introduced by Newman and Girvan|71 123| 
and are also similar to the networks studied by Watts at the end of last 
century . A particular case of this networks is the one shown in figure . 

In a random network, clustering is equal to the probability that two agents 
are acquainted, motivating us to define guild clustering as the clustering co- 
efficient of a guild in the absence of the rest of the network. We studied the 
dependence of influence on guild clustering (fig. l|HP)) finding that there is a 
linear dependence with a slope that decreases with a, meaning that the rein- 
forcement mechanism is more effective in sparsely connected networks than in 
strongly connected ones. Either way, the effect continues to remain mild com- 
pared to its scale-free counterpart. On the other hand, the characteristic time 
T required for the system to converge increases with decreasing aasT~ a~ K ^ b ' 
with dn/db < (see Fig. Ef>). 



5 Equilibrium Distribution 

Using the formalism recently introduced we can also find how states are dis- 
tributed when the network reaches equilibria. In order to do this, first we notice 
that eqn.Q vanishes when 

T l pi 

If we multiply eqn. (0J by Pji and add over all i we will find a rate eqn. for 
Yj Q . This, must vanish in equilibria implying that 

^E^f = E T ^ ( 12 ) 

i i 

The term on the righthand side does not depend on i so the summation just 
yields the number of guilds we have chosen to approximate the system. Con- 
sidering a system with m guilds leads us to conclude that (|T2|) is the same 
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as 



1 



E 




(13) 



m 



The lefthand side of this equality does not depend on j, so when the system 
reaches the steady state, the ratio between T Jct /Tj must be the same for every 
guild. Combining this with condition (|11[> we conclude that upon equilibria the 
system reaches a configuration in which the ratio between the number of agents 
in a certain state over the total number of agents is independent of the degree, 
thus each guild has exactly the same distribution of states as the entire network. 



6 Network Structure 

The formalism presented above can be used to understand some topological 
characteristics of the studied network. Recently, two groups have presented 
works indicating that a simple dynamical process like the one studied in this 
paper can be solved, in the case in which two states are involved, by the use of 
conserved quantities or martingales |25l |2E] . The first of them proposes that a 
weighted average over degree should be conserved in a dynamical process of this 
kind. In our formalism we average the fraction of agents in each state weighted 
by their relative degree as 



We can shows that this is a conserved quantity by taking eqn. (J2J, multiply it 
by Pji and add it over i, j and a to get, 



Performing the summation over a on the first term and over /z in the second 
term and cancelling the Ti we get. 



which is clearly equal to zero. This quantity is naturally conserved if we consider 
that Ti a represents the fraction of agents in the a state that are known by 
someone on the the i'th guild. Thus Ti represents the total numbers of agents 
know by the i'th guild, and T therefore represents the fraction of agents known 
by all the guilds combined, or in the other words the number of non-intersecting 
sub-graphs of the network that are able to reach all the nodes of it. 
In a two state system in which each one of them is represented by a spin up 
or down we can show that the degree averaged magnetization does not depend 
on time by replacing the fraction of agents in a certain state by the fraction of 
agents in that state times its spin. 




(14) 





Pia. 
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If we span Ti a we can show that the change on this quantity can be constructed 
by using our formalism as 



where the two added terms are strictly equal if we consider that [i and a are 
mute variables. This agrees with the results found by Wu et. al. [55] and 
can be extended to an arbitrary number of spins. Which in this case are noth- 
ing more than labels, and we believe its names should not enter the calculations. 



7 Conclusion 

Using a topologically considerate mean-field approach we were able to charac- 
terize the simple dynamics of a node updating mechanism on a complex network 
with an arbitrary degree distribution. 

The main idea behind this technique was to group nodes with similar topological 
characteristics into aggregates we called guilds. We characterized the reproduc- 
tive capability of guilds by a quantity we called influence and showed that in 
the case of a scale-free network the system bias dramatically towards the states 
pertaining to the highest degree guilds as the characteristic exponent decrease, 
while on an exponential network this effect is comparatively mild. 
It was also shown that the convergence time scale as a function of connectivity, 
with a characteristic exponent close to 1, when we considered a network formed 
by mildly connected clusters. This is a consequence of the reinforcement mech- 
anism introduced by the node updating rule. We were also able to show that 
clustered groups tend to be more influential than un-clustered ones. This is due 
to the fact that they are harder to invade, allowing them to act over the system 
for a longer period of time. 

The final configuration of the system was unravelled, showing that it arrives to 
a stationary distribution in which the fraction of nodes with a certain degree 
that are in a certain state, over the total number of nodes with that degree, is 
independent of the particular degree and it therefore mimics the network. 
In the last section we discussed a purely topological quantity, which is natu- 
rally conserved and represents the number of subgraphs present in the network 
with the ability to held it together. This suggests a direct link with network 
resilience and will be present in future discussions. We also showed that the 
degree weighted average of states over the system remains constant, regardless 
of the value and number of them in agreement with previous results. 

Work supported in part by Fondecyt project 1020829. PAM acknowledges 
the support of the Santa Fe Institute through and International Fellowship. 
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